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Abstract It is argued that, in order for the gravitational field to be 
propagated as a wave, it is necessary for it to satisfy a further set of field 
^^"^ equations, in addition to those of Einstein and Hilbert, and these equations 

i mean there is a preferred coordinate frame, called the Global Inertial Frame, 

I giving rise to a unique metric . The implication is that a true gravitational 

' field is not compatible with Einstein's Principle of Equivalence, which is in con- 

tradiction with his other fundamental concept of locality. The additional field 
equations ensure that gravitational collapse does not go below the Schwarzschild 
radius, thereby excluding the possibility of singular solutions (black holes) of the 
Einstein-Hilbert equations. Such solutions would also violate Einstein's locality 
principle. 
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u ' 1 Introduction 

In 1907, in two articles[l][2] separated by a few months, Albert Einstein stated 
two principles to which he would subsequently give the names Principle of Local 
Action (PLA), or in the original Prinzip der Nahewirkung, and Principle of 
Equivalence (PE). The first of these states that nothing goes faster than light, 
and it was considered by him to be the principle underlying what we now call 
Special Relativity, Einstein's theory of 1905. He restated it many times during 
his later life[3][4][5], especially in his criticism of quantum theory, and in the last 
of these citations he stated the view that it would be impossible to do science 
' without it. As for PE, although he classified it as "the happiest thought of my 

, life" ([6] pl78), leading him, as it did, to the 1915 theory which we now know 

as General Relativity (GR), the statement of 1907 was very modest compared 
with what it later became; it stated that no observer can distinguish between a 
uniformly accelerated frame and a frame at rest in a uniform gravitational field. 

Even before the creation of GR, the two principles presented some com- 
patibility problems; both the gravitational red shift and the variation of the 
free-space refractive index, predicted by Einstein in 1911 and resulting in the 
bending of light beams, made the limiting value of c somewhat ambiguous. But 
the search from then to 1915, by Einstein in collaboration with Grossman[6] and 
independently by Hilbert [7], for a group of transformations more general than 
that of Lorentz under which the laws of field transmission should be invariant, 
led eventually to a very strong formulation of PE. From that point on Einstein 
was to state repeatedly that all the laws of physics should be independent of the 
coordinate system. This strong form of PE is, I shall argue, totally destructive 
of PLA. 

During the formative years of GR he and Grossman tried to create a theory 
of gravitation based on a more restricted group of transformations [6] than is 
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demanded by PE. An examination of the articles they wrote during this period 
shows that it was a desire to incorporate PLA, in the form of a tensor represent- 
ing the local energy density of the gravitational field, that motivated the various 
twists and turns taken by them. Einstein had decided by 1915 that this was 
not possible, and in articles by Hilbert[8] and Schr6dinger[9] the impossibility 
of constructing such a tensor was confirmed. Nevertheless, he returned in 1918 
to the idea that gravitation, like any other field, must propagate at the velocity 
c. By imposing on the gravitational potential, which in conformity with PE is 
the same as the metric of the curved space, a certain restriction, he deduced a 
formula[10] for the total radiation emitted even though the local flux thereof re- 
mained ambiguous. It was not long before Eddington[ll] drew attention to the 
fact that Einstein's waves could be transformed out of existence by a coordinate 
change permitted imdcr PE; he observed caustically that the waves "travel at 
the speed of thought". However, both Hilbert and Eddington ([11] pp40-41) 
showed themselves ready to modify PE in order to accommodate some notion 
of locality. The former [7] actually identified the related notion of "causality" 
as a further requirement for a physical theory of gravitation, and he stated a 
criterion to meet this requirement, that the signature of the metric gij must be 

preserved as (H ) throughout any evolution of a system. I shall show that 

Hubert's causality is a close relation of Einstein's PLA, and that playing fast 
and loose with the signature of gij is precisely what has led a substantial part 
of our community into the blind alley of black-hole theory. It should be noted 
that more recently[12] it was shown that coordinate transformations exist which 
transform away the Einstein energy loss globally as well as locally. A strong PE 
is not compatible with the existence of gravitational waves. 

With respect to black holes Oppenheimer and Snyder[13] showed in 1939 
that a fairly simple metric describes what may plausibly be considered a set 
of dust particles under zero pressure which, in a finite time interval, collapses 
into a black hole, so called because light cannot escape from it. Objects with 
nonzero pressure are now widely believed to undergo similar collapse, subject to 
certain conditions of size and internal temperature, but the dust cloud remains 
the basis of the black-hole paradigm. It seems likely that Einstein knew of this 
work before it was published, because almost simultaneously he published an 
articlc[14] criticizing the concept, on the groimds that the crossing of the event 
horizon by the dust particles would violate his locality principle. In this he was 
also affirming his support for Hilbert's notion of "causality", referred to above; 
at the point of crossing the event horizon, the signature of the metric changes 

from (H ) to ( — I ). In other words PE, a principle which Einstein 

apparently still supported, predicted a process which violated locality; it seems 
that locality won the day as far as he was concerned! 
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2 The relativistic theory of gravitation 



GR has only one set of field equations 

n'^ -\g''gkin''' = iT^KT' , « = ^ , (1) 

where TV^ is the contracted curvature tensor derived from the metric tensor gij , 
and T^^ is the local material stress, or energy density tensor. In his article on 
gravitational waves[10] Einstein introduced the noncovariant condition 

di<P'^=0 , $^^=5^^V^ . (2) 

This is often referred to nowadays as a gauge condition, and it was further 
developed by de Bonder [15] and later by Fock[16], who showed that, if the 
coordinates are cartesian and satisiy (2), then they also satisfy the harmonic 
condition 

ax, = o , (3) 

where □ is the generalized d'Alembertian operator 

□ = ^5,y'^/^afc . (4) 

Fock showed that many complex calculations in GR are greatly simplified in 
harmonic coordinates. In particular he used them to calculate the rate at which 
a binary system, of which the only example for which he then had the relevant 
data was Sun-Jupiter, radiated away energy through gravitational waves. Al- 
though the effect was too small to observe in that system, his calculation formed 
the basis for such a system involving a neutron star, which made possible the 
calculation [17] of the radiative energy loss in a binary pulsar. 

Although Fock criticized PE, and indeed considered that it was incorrect 
to label the Einstein- Hilbert theory " General Relativity" , he did not consider 
the harmonic system of coordinates as in any way privileged. That was left to 
a later school of gravity theorists in the Soviet Union[12] who proposed that, 
within a given family of metrics classified as equivalent under PE, it is the har- 
monic one, which correctly describes the physical system. There is still a family 
of coordinate franicjs with this metric, but the group of coordinate transfor- 
mations connecting these frames is the Lorentz group of the Special Theory of 
Relativity. They called the resulting theoretical structure the Relativistic The- 
ory of Gravitation (RTG). A gravitating system studied in isolation, for which 
the distant field is Newtonian, is classified as an island system, and I shall call 
the harmonic frame for such a system, in which its centre of mass is at rest, the 
global inertial frame (GIF). 

RTG differs significantly from General Relativity. In particular gravita- 
tional collapse does not go below the Schwarzschild radius, and so there are 
no black-hole singularities. At the same time the theory includes unambiguous 
expressions for the material and gravitational stress tensors, and the latter is 
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associated with a real physical field which propagates at the speed of light, just 
like the electromagnetic field of Faraday and Maxwell. The quantities are 
considered to be gravitational fields instead of being the gravitational potentials 
of GR, and the field equations of RTG consist of (1) together with (2). It is 
a radical departure to treat this latter equation as an essential, globally valid 
field equation, which, for a given source, results in a unique field, defined with 
respect to the GIF. 

A simple example of an island frame arises from considering the Schwarzschild 
metric for a point mass m at the origin 

^^2^r-2m^^2 r ,^2_^2r^g2^^.^2^^^2^^ 

r r — 2m ^ ' 



for which the d'Alembertian is 
r — 2m r 



□ = _^ df - ^ [r (r - 2m) + 2{r - m)dr + dg + coiOde + csc^ 6»5|] . 



(6) 

The coordinate z = r cos is not harmonic, while Z = (r — m) cos 6 is, so the 
metric associated with the GIF is 

ds^ = ^^de - ^^^dR^ -{R + mf (de^ + sin2 edcj>^) . (7) 
R-\- m R — m 

Using this metric instead of Schwarzschild's gives a small correction[12] to the 
perihelion advance of a planet's elliptic orbit, to the bending of light, and to the 
gravitational red shift, but they are all far too small to be measured. 

3 The gravitational collapse of a dust cloud 

Oppenheimer and Snyder[13] studied the metric 

ds^ = dr^ - S\dR^ - {dO^ + sin^ 6»#2) , (8) 



where 



and 



Note that 
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(9) 



1 (i? > 1) 



yi?75 [R > 1) 

S/R (R < 1) 



(11) 



which means that this metric has a discontinuity at i? = 1. They showed that 
a frame defined by it is comoving in the sense that its material stress tensor, 

T^' = ^ U'' - Ig'^gkiTZ'^') , (12) 
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has only a timelike nonzero component, that is 



rpcxf3 rpaO r\ 'TiOO 

for which the mass density is 

pdRdedcj) = T°° yf^dRdedcj) = ^ sin OdRdedcj) . (14) 

27r 

They used this metric to model the gravitational field of a spherically symmetric 

dust cloud with no rotation; the function / (i?) contains the i?-dependence of p, 
and in the exterior region, i? > 1, it is constant, making p zero. The comoving 
property of this frame means that the timelike coordinate r is a local proper 
time. It is a simple matter to establish that each particle in the cloud reaches the 

"event horizon" , a point at which the (H ) signature of space-time changes, 

within a finite interval of r. This was taken in [13] to show that the cloud 
collapses into a black hole, but we shall show, using the harmonic coordinates, 
that this was an incorrect conclusion. 

For the exterior region, i? > 1, the latter authors looked for a coordinate 
transformation (r, R) — > {t, r) such that the above metric transforms into the 
Schwarzschild metric with 2m = 1 . This leads to the result 

f=^(i?3/2_53/2j +2-2^5+2 log -log (S-1) , S = r .(15) 

They then obtained an interior metric by imposing the conditions g*^ = and 
S = r there also. They required in addition that t, but not Ir, be continuous 
at i? = 1, and they deduced that 

i = -log(t/-l) + 21og^^ + ^-|y3/2_2^ ^ (16) 

where 

R'^-l r 

y=^^R ■ ^''^ 
The coordinates (t, r) they refer to as " external coordinates" for both the exte- 
rior and the interior regions; this is by contrast with the comoving coordinates 
(r, R) . In the external frame t becomes infinite for the exterior at r = r^o = 1 
and for the interior at 

roo = lR-\R^ {R<1) , (18) 

which is the event horizon. 

All the particles of the cloud go to their final positions, given by Too (R) and 
stay there; this is a stationary state and it is surprising that Oppenheimer and 
Sneider did not acknowledge it as such. What happens on the wrong side of the 
event horizon, that is for proper times satisfying 

2 2/3-i?2\^/^ 



5 



has no physical relevance at all. What seems to have caused their confusion was 
the discovery that the metric coefficients become singular near r = 1. There is a 
simple explanation for this, since the Schwarzschild metric in the exterior region 
has the same behaviour. Furthermore, the density p{r co) becomes infinite at 
Too = 1 , as may be shown from the constant cumulative density in the comoving 
frame, namely 

P{R)= A-Kp{R)dR = R^ . (20) 

Jo 

The corresponding quantity in the external frame is 

P (roo) = [R (roo)]' = 8 cos3 [- + , (21) 

and near Too — 1 this gives ciii infinite density, namely 



47rPooW~\/^^ • (22) 

The infinite value of Qrr found by Oppenheimer and Sneider in the limit r ^ 1 
is simply a consequence of this mass concentration at the surface of the cloud. 
By contrast gtt is zero at all points on the event horizon, which indicates an 
infinite red shift in the limit t ^ oo. 

Their exterior solution (15) is not unique, owing to the arbitrary function in 
their eqn.19 having been equated to one. Of course the absence of uniqueness is 
inevitable in General Relativity, it being a manifestation of PE, which ordains 
that no coordinate system is to be favoured. In RTG, by contrast, the external 
solution is unique, because it has to satisfy the harmonic condition; in that case 
the external solution is identical with (15), except (see the previous section) 
that S is equal to r + 1/2 instead of to r. As for the internal solution, the 
condition g*^ = of Oppenheimer and Snyder somewhat fortuitously ensures 
that tji/rji is continuous at i? = 1 if r and t are, and this then makes gu and grr 
continuous there, even though guR is not. The condition g*'^ = is arbitrary 
though allowable under PE, but it is not satisfied by the RTG solution. That is 
constructed by imposing the harmonic condition, and by requiring the continuity 
of both tr and tji at R = 1, thereby setting a constraint on the function f{R) 
in the internal region; this ensures by a different route the continuity of the 
metric. 

In accordance with the latter condition, I choose for the interior density 
function 

/ {R) = R^/^e^^/^ , X = l-R {R<1) (23) 

giving 

Sr = S^ , + (i?<l) , (24) 

and the mass density 

ovp2^3X 

pdRdOd(t>= sin6dRded(t> . (25) 



6 



The d'Alembertian operator in terms of the coordinates {R, S, 6, (j)) is 

1 



□ = Di - ^ [dl + cot ede + csc2 edl) , (26) 



where 



1 



ds 



2 / r 



S_ 

Or . (27) 



The harmonic coordinates 

ajj = {t, r sin ^ cos ^, r sin ^ sin ^, r cos ^) 

satisfy 

□it = , (28) 

and 

°i + ;^)^ = o ' (29) 



for which the exterior solutions, obtained by putting / = 1 and ^ = y/R/S^, 
are the same as (15) for t, with r replaced by 

r = S-^ . (30) 

This is, as expected, the harmonic form of the Schwarzschild solution, as ob- 
tained already in the previous section. The interior coordinates are obtained 
by solving these same equations with r,t,rR and Ir continuous at i? = 1, and 
t and r having the Correspondence Principle (CP) behaviour t ^ T,r ^ S as S 
tends to +00. The latter conditions give the leading term for large positive S, 
namely 

f^_253/2g3Z/2^_2^3/2g3Z/2 £"^ = 6^7? , (31) 

3 3 . V ; 

or equivalently 



it 



2/3 



r~roe-^ , ro = (^-^j , (32) 
leading to a density 

3^2 

p (r) sin edrdedcj) - - — ^ sin edrdOdcj) {r < ro) • (33) 

47r'rg 

This gives us the collapse of a sphere of radius rg {t) with initially uniform mass 
density, and we see that the uniform density is preserved during the first stage 
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of collapse, in agreement with the Newtonian behaviour for such a density. The 
leading relativistic correction is obtained from the next term, namely 



t ^ _253/2g3Z/2_^l/2 ^ ie-3Z/2\ 

3 V 2 2 / 

from which S may be eliminated to give 



S--e-^ + -e-^^ (34) 
4 4 



,.e-(ro + i-le--) , = (-|) - ^ , (35) 



p (r) sin edrd0d(j) - ^^-3 1 - h ) sin Odrdedcj) (r < ro) . (36) 



with the density 

3r2 A_ 3 5r 
47rro V 4ro 4rQ 

This latter result shows a departure from the uniform density associated with 
Newtonian gravity; as the collapse progresses the density near the surface (r = ro) 
increases faster than at the centre of the cloud. 

To follow the evolution into the region of strong gravity requires a numerical 
approach, for which we make the further change of variables from R to X and 
S toY = \ogS + Z, leading to 



5^ 3 3 2 + 3X 



dx + ^dy ] (37) 



From this it is possible to derive more terms in the asymptotic series for r and 
t, namely 



(38) 



and 



r~e-^(^e^+X:«ne-''^/2j , (39) 
where the first few nonzero terms are 
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2 3X/2 

3 

2 7 15 1 

-e^ H e 

5 8 4 



-2Z 



-4Z 



40 



ie-3^/2 + ^e^/2(l3-6X + X2) 



20 



43 



29 



— e H e~ e" 

70 96 10 32 



-22 + + ^( 

160 3360 



-6Z 



(40) 
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and 



■"3 



-7 + 7^"'^ 
4 4 

2 2 ^ 



6 4 ^ 



4X + 9) 



_^ -3X , 100 3Ji:/2 

12 3 



(41) 



The partial differential equations (28) and (29) must now be integrated sub- 
ject to the boundary conditions at X = 0, for which the first equation displays 
a logarithmic singularity at 5 = 1, that is y = 0. This suggests there is a 
similar singularity along the characteristic Y = Yq {X) through (0,0). This is 
in fact the event horizon in the {X, Y) coordinates, and it is determined by the 
ordinary difii'erential equation 

Y^ = iJ^e^->'o/2 _ g-3X/2-3Yo/2 

= Xe^-W2 + (i_J^)e-X/2-2Fo_g-3X/2-3ro/2 ^ (43) 

which integrates numerically to give the function 



X 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


Yo 


.005 


.021 


.049 


.091 


.150 


.227 


.326 


.447 


.591 


.758 



It may be noted that this is finite at X = 1, even though some coefficients 
of the partial differential equations are infinite there. This is a consequence of 
the choice made for the variable Y. A further advantage of this choice is that, 
as y — > +00. the characteristics are Y =constant, which aids the numerical 
integration. The latter is achieved by replacing each partial differential equation 
by a set of ordinary differential equations for the vectors r (X) = r {X, Yi) and 
t (X) = t {X, Yi) , where Yi are a set of 100 values of Y spaced equally between 
an upper value Yq at which the above asymptotic expansions apply (I took 
lo = 5) and a lower value Yi close to zero (I took Yi = .001). The integration 
confirms that t becomes infinite as Y approaches lo (X), which means that we 
are seeing an extreme form of the gravitational red shift. The values of r {X) 
along this same curve, which we designate Too (X), are 



X 





0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


^00 


.500 


.500 


.500 


.499 


.497 


.492 


.481 


.456 


.398 


.272 






and this shows that gravitational collapse ceases at the radius rooo, which 
takes the value 0.5 in our units. The latter fimction, or rather its inverse X (roo), 
gives us the limiting density (r), since the density in terms of the comoving 
coordinate X remains constant throughout, so that the cumulative density is 

Poo (r) = r 47rp^ (r') dr' = [1 - X (r^)]' e^^^''-) . (43) 
Jo 

A corresponding uniform density would give 

Puir) = ^ = 8rl , (44) 

' ooO 
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and a comparison may be made by plotting Poo/Pu a-s a function of r. This 
shows that the process indicated by (36) above, namely a concentration of par- 
ticles at the surface, continues and intensifies. For example, we find that Pu 
would give less than 2% of the particle distribution in the range .497 < r < .5, 
as compared with Pqo giving more than 28%, while, in .481 < r < .5, Pu gives 
less than 11%, compared with Poo giving more than 61%. As Too tends to 0.5 
the derivative of Poo, that is becomes infinite; this is natural, because the 
metric coefficient Qrr must become infinite at r = 0.5 to match with its value 
in the exterior region P > 1. Also, associated with the infinite red shift at the 
event horizon, the determinant g of the metric approaches zero there. 

4 Conclusion 

Our result shows that, contrary to the general result claimed by [13], there is, 
for a suitable choice of the density function / (P) , and of course a suitable 
coordinate frame, a stationary solution of the Einstcin-Hilbcrt equation, and 
this state is approached in the limit t +00. It should be noted that, both 
in the present article and in [13], an additional condition had to be imposed in 
order to obtain an interior metric satisfying the Correspondence Principle, and 
that in both cases there is a singularity in the metric in the limit t +00. 

More significant, however, is that the insistence on a global inertial frame as 
the vehicle for gravitational waves restores the gravitational theory of Einstein 
to the framework of his Special Theory. In particular the waves travel with the 
velocity of light and are consistent with his principle of locality. It is entirely 
natural then that RTG should also rule out gravitational collapse beyond the 
event horizon, leading, as it would, to an exotic topology in which the space and 
time coordinates would change places and to a complete demolition of causality. 

The key to understanding what happens as a particle approaches the event 
horizon is that such an intense gravitational field produces a red shift which 
becomes infinite, that is all physical processes including gravitational collapse 
are infinitely slowed down. Effectively time is frozen. There is also an infinite 
mass density at the surface of the cloud. This is an extreme situation which 
none of us is likely to experience directly, but it ensures that, even in such 
circumstances, the world remains local and comprehensible. 
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